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Abstract 
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lated.  The  use  of  a  laser  to  enhance  surface  charge  density  is 
discussed. 
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1.  Introduction 


Recently,  there  has  been  much  interest  in  the  effects  of  laser 
radiation  in  solids.^"  Using  experimental  electronic  band  struc¬ 
tures  and  time  dependent  perturbation  theory,  researchers  have 
calculated  the  photon  absorption  rate  for  the  electronic  states 
in  a  variety  of  semiconductors.  Much  effort  has  also  been  devoted 

to  the  study  of  the  effects  of  laser  radiation  on  the  phonons  in 
2 

solid  surfaces.  However,  the  effect  of  a  surface  on  photon-excited 
electrons  has  not  been  examined. 

Studies^  of  the  effects  of  synchrotron  radiation  on  metal  sur¬ 
faces  have  shown  that  photon-stimulated  desorption  can  occur  through 
electronic  excitation.  The  details4  of  this  desorption  process  are 
involved,  but  the  ultimate  reason  for  the  desorption  is  the  shift 
of  electronic  charge  and  the  resultant  Coulombic  repulsion. 

For  a  semiconductor,  the  bulk  electronic  band  structure  con¬ 
sists  of  a  number  of  valence  and  conduction  bands.  In  addition 
to  these  bulk  bands,  a  surface  can  introduce  additional  bands  and 
local  states.^  These  surface  bands  correspond  to  a  charge  density 
localized  in  the  surface  region.  Consequently,  by  using  a  laser 
to  excite  electrons  to  or  from  the  surface  band,  we  can  effectively 
control  the  surface  charge  and,  ultimately,  the  Coulombic  repulsion 
or  attraction  of  the  surface  for  adsorbed  species. 

In  the  following  section  we  will  examine  the  effects  of  a 
surface  on  the  band  structure  of  a  model  system.  Laser  excitation 
of  valence  electrons  to  surface  states  in  this  model  system  will 
then  be  analyzed  in  the  next  two  sections.  Finally,  a  discussion 
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of  the  applications  and  limitations  of  this  model  for  studying 
laser  induced  surface  processes  will  be  presented. 

2 .  Model  Band  Structure 

In  an  infinite  crystalline  solid,  the  energies  of  the  elec¬ 
trons  divide  into  bands  of  allowed  states  separated  by  energy 
gaps.  For  a  finite  crystal  thin  film  calculations®  have  demon¬ 
strated  that  energy  states  exist  in  these  gaps.  Furthermore, 
the  charge  associated  with  these  states  resides  mostly  in  the 

surface  region. 

7  5 

Heine  and,  later,  Lundqvist  have  shown  that  truncating 
a  linear  chain  produces  similar  effects  on  the  energy  gap.  For 
this  simple  one-dimensional  model,  the  dispersion  relationship 
is  given  in  Figure  1.  In  a  semiconductor,  the  lower  inverted-U 
curve  would  be  the  bulk  valence  band  and  the  upper  U-shaped  curve 
would  be  the  bulk  conduction  band.  The  curve  in  complex  crystal 
momentum  space  is  associated  with  the  surface  states.  If  we  have 
a  chain  with  lattice  constant  a  parallel  to  the  z  axis  in  the 
region  z<  a/2,  the  interior  wave  function  associated  with  these 
gap  states  is  that  of  a  damped  oscillator: 

ipK(z)  *  Cse~<z  cos(^  +  <f>)  (1) 

with  corresponding  energy 

E{<)  -  (|)2  —  k2  ±  (Vg-g2K2)1/2  (2) 

where  C  is  the  normalization  constant;  k  is  the  imaginary  part 
s 

of  the  crystal  momentum  and  ranges  from  0  at  the  edges  of  the 
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gap  to  -|Vg/g|  near  the  center;  g  is  the  reciprocal  lattice 
constant  and  equals  2ir/a;  and  V  is  the  g  component  of  the  Fourier 
transform  of  the  effective  potential.  The  phase  factor  $  is  given 
by 

sin  2<J>  =  $S  .  (3) 

g 

The  interior  bulk  states  would  have  the  form 


4>k(z) 


C[eikz 


+ 


E (k) -k2  i (k-g) z, 
V  J 

g 


(4) 


with  corresponding  energy 


E (k)  =  |[k2  +  (k-g)2  ±  (  [k2-(k-g)2]2  +  4|v  j2)1/2]  (5) 
z  y  • 


The  interior  solution  for  both  the  surface  states,  equation  (1) , 
and  the  bulk  states,  equation  (4),  should  be  matched  to  the 
exterior  solution  of  the  form 


ip(z)  =  C  e'qz 

S 


(6) 


with 


q  =  /Vq-E(k)  ,  (7) 

where  VQ  is  the  work  function  plus  the  energy  at  the  top  of  the 

valence  band.  For  bulk  states,  C  in  equation  (6)  and  E(k)  in 

s 

equation  (7)  should  be  replaced  by  C  and  E(k),  respectively. 

To  obtain  the  normalization  constants,  we  assume  the  chain 
to  be  of  length  L =  Na  and  the  charge  density  was  symmetric  about 
the  center.  Under  these  conditions  we  obtain 


C 


{[1  + 


.E  (k)  -k4. 
V  V  ' 


]  L  + 


1 

7 


(8) 
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for  the  bulk  states  and 


C  -  te-«a(2*  COS  2f-asin  2j>  (1  + 
S  (2c)  1  +  q 


,  kL  ,  , 

+  le-  Sll  + 

2k 


,-qa  .  I 
T“]  2 


for  the  surface  states. 

In  the  following  sections,  we  wish  to  examine  laser  excita¬ 
tion  from  the  bulk  valence  band,  equation  (4) ,  to  the  surface 
band,  equation  (1) .  The  initial  conditions  are  assumed  to  be  a 
filled  valence  band,  but  empty  surface  and  conduction  bands. 
Since  equation  (4)  represents  electrons  delocalized  throughout 
the  chain  and  equation  (1)  corresponds  to  electrons  in  the  sur¬ 
face  region,  the  excitation  of  electrons  to  the  surface  states 
will  increase  the  surface  charge. 

3.  Selection  Rule 

To  determine  the  transition  rate,  we  must  evaluate  an  in¬ 
tegral  of  the  form 


M  =  <K|A*p|k>, 


where  A  is  the  vector  potential  of  the  laser  radiation  and  p 
is  the  momentum  operator  of  the  electron.  Under  the  dipole 
approximation  and  assuming  the  laser  is  polarized  parallel  to 
the  chain,  M  becomes 

OirT  1  /  ^  j 

«  ■  -  ilm’  «l§llk>'  (11) 


where  I  is  the  laser  intensity  and  oj  is  the  angular  frequency. 
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The  wave  functions  of  interest,  equations  (1)  and  (4) ,  can 
be  written  in  the  form 


*k(z)  =  e_lkzuk(z)  (12) 

<PK(z)  =  e~x{2)zeKZ  uK(z),  (13) 

where  the  functions  uk(z)  and  u^z)  have  the  periodicity  of  the 
lattice.  The  integral  in  equation  (8)  now  becomes 

M'  =  «||j|k>  =  2/dz2  e"l(k~2)  zeKZ[uK(2)u]^(z)-ikuK(z)uk(z)  ] 

21  L 
2  -  2 


l~°  "f(l+2fc) 


[^(zju^fzj-iku^fzjv^tz)  ]  ,  (14) 


where  we  'have  assumed  that  for  large  L  the  contribution  of  the 
exponential  tail,  equation  (6),  to  the  integral  will  be  negligible. 
If  we  change  the  variable  of  integration  to  z1  =  z  +  Ha  and  exploit 
the  periodicity  of  uk(z)  and  u^(z),  we  obtain 

M*  -  NI1e1(k~l)  Za  eK^a/dz^e-1  (k~^}  z' eKZ‘  [u  (z 1 )  u/  ( z 1 ) 

«,=0  a 

2 

-iku^tz’Ju^tz’)  ] 


V  ei(k-f)ia  jcla  ,  ±-lk> 

do  e  2  e  '  dz'  k  o' 


(15) 


where  the  subscript  "o"  indicates  integration  over  the  first  unit 
cell.  The  transition  probability  is  related  to  the  square  of 


equation  (15)  : 


M'2  =  I  Z  ei  ae<(£+m)a  <K  ||_jk>  .  (1 

m=0  =0  dz  ° 

If  we  change  the  range  of  summation,  the  sums  in  front  of  the 
integral  become 


N-l 

2 


S  =  eKNa  Z 


N-l 

2 


£  ei  (k-J)  (2,-m)  ag<  (i +m)  a 


(17) 


N-l  5_  N-l 


Changing  the  summation  variable  to  s  =  £  +  m  and  d  =  £-m,  we  get 


S  -  V  ei(lt-?lia  "1^  e*sa  (18 

d=l-N  s=l-N+d 

where  the  prime  signifies  a  count  with  an  increment  of  2.  If 
N  is  very  large,  the  sum  over  s  will  only  be  significantly 
affected  by  large  values  of  d.  Therefore,  for  large  N,  we  can 
approximate  this  function  by 


S  *  e 


icNa 


N-l 


Z  e 
d— (N-l) 


,i  <N;)> 


Z  ’  e 
s=- (N-l) 


Ksa 


(19) 


Taking  the  limit  as  N  goes  to  infinity,  the  function  takes  the 
form 


S  z  J^6(k-f) 
a  2 


1-e 


2k  a 


(20) 


where  <5  (k-^)  is  the  Dirac  delta  function.  Equation  (16)  can 
now  be  written 


M’2  s  6(k-|) 


1-e 


2<a 


(21) 


Consequently,  the  transition  from  a  bulk  to  a  surface  state  is 
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only  permitted  if  the  real  part  of  the  crystal  momentum  remains 

unchanged.  This  relationship  is  not  too  surprising  since  it  is 

an  exact  restriction  on  laser-induced  transitions  between  bulk 
8 

bands.  Furthermore,  for  our  model,  it  confines  us  to  the  top 
of  the  valence  band,  where  the  density  of  states  is  a  maximum 
(infinite)  and  the  laser  frequency  needed  for  a  transition  is 
a  minimum. 


4 .  Transi tions 

To  first  order,  the  transition  rate  from  a  valence  band  to 
a  surface  band  is  v 

.  g 

T  =  2uh  (|2)  /dk/dK  |  <k  |  A*p  |k>  |  26  (E(k)  -E  (k)  -<i>)  .  (22) 

g  0  0 

Using  equations  (11)  and  (21) ,  this  can  be  simplified  to 

V 


„  2  — *  I  <<  |~|  2.?  I2 

T=  (TTf)  (^SLy)  JdK?  - dZ  La  S  (E  (<) -E  (k) -u>)  . 

V  ai2  0  1  -  e^ 

q 


(23) 


After  evaluating  the  integral  over  k,  the  expression  reduces  to 

d<  i 


137'  '  ,  2'  .  2<a  1  dE  (<) 

v  iii  i  —  e 


(24) 


where  k  in  equation  (24)  refers  to  the  state  obeying  the  resonant 
condition: 


a)  =  V  -  k2  ±  (V2  -  g2<2)1/'2. 

g  g  * 


(25) 


We  now  define  the  absorption  cross-section,  o: 


_  _  wT 

a  =  T  ‘ 


(26) 
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Using  equations  (1),  (4),  (8)  and  (9)  in  equation  (24)  and  tak¬ 
ing  the  limit  as  L  goes  to  infinity,  the  cross-section  becomes 


With  equation  (2),  we  can  also  readily  evaluate  the  derivative: 


Equations  (27)  and  (28)  constitute  the  cross-section  for  electronic 
transitions  from  the  valence  band  to  the  surface  band.  Although 
this  cross-section  is  quite  complicated,  we  can  readily  deduce 
its  behavior  by  analyzing  the  expressions  at  various  limits. 

If  the  exciting  laser  radiation  is  at  the  frequency  near 
^  Eg  where  k  =  -|Vg/g|,  equation  (28)  will  vanish,  and  thus 


0M. y  •  0.  (29) 

*  g 

This  is  exactly  what  one  would  expect  since  this  mid-gap  energy 
is  a  branch  point  at  which  no  surface  states  exist. 


If  the  laser  radiation  is  near  a  frequency  of  0  or  Eg ,  the 
cross-section  becomes 


2  E  +g  9  , 
9  Eg-9 


and 

0<i)-*-E 

g 


E  [E  +g2  ] 

g  g  y  1 


(31) 
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At  both  extremes,  k  -*■  0  and  equations  (30)  and  (31)  diverge. 
This  occurs  because  at  the  surface  band  edges  the  charge 
associated  with  the  surface  states  becomes  more  and  more  de¬ 
localized  throughout  the  lattice  until  at  k  =  0 ,  the  charge 
is  completely  delocalized.  At  this  point  the  surface  states 
become  bulk  states,  and  instead  of  cross-sections,  one  should 
consider  absorption  coefficients. 

Figure  2  illustrates  the  behavior  of  the  cross-section 

g 

over  the  entire  frequency  range.  The  values  for  the  lattice 
constant,  a,  and  the  energy  gap,  E^,  were  taken  to  be  those 
of  silicon. 

5 .  Discussion 

As  pointed  out  in  the  Introduction,  we  are  interested  in 
using  a  laser  to  increase  the  surface  charge  density.  Since 
the  effective  charge  depth  is  |l/(2tc)|  (see  equation  (1)),  we 
wish  to  excite  states  for  large  values  of  < .  For  silicon  at 
maximum  k,  the  surface  charge  depth  is  1.39  lattice  constants 
(1.39  a);  consequently,  the  charge  is  confined  to  a  region 
very  near  the  surface.  As  we  move  away  from  this  mid-gap 
region  the  charge  depth  increases  as  mentioned  in  the  last 
section.  Therefore  we  are  mostly  interested  in  laser  frequen¬ 
cies  near  ^Eg-  As  seen  in  figure  2,  the  cross-section  is 
quite  substantial  near  the  mid-gap  region  and,  subsequently , 
we  could  readily  increase  the  surface  charge  by  using  a  laser 
of  moderate  intensity.  Even  at  to  =  .15  E  and  w  =  .85  E  ,  where 
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the  cross-section  has  increased  substantially ,  the  charge  depth 
is  still  only  about  2  lattice  constants.  Consequently,  we  would 
expect  a  laser  tuned  to  a  frequency  around  j E to  be  an  effective 
controller  of  surface  charge.  Furthermore,  since  the  energy 
gaps  for  most  semiconductors  are  approximately  1  eV  (for  silicon, 
E^=  1.17  eV)  ,  the  frequency  of  interest  will  be  in  the  infra¬ 
red  region. 

Since  these  results  are  based  on  a  one-dimensional  model, 
certain  limitations  should  be  pointed  out.  First  of  all,  most 
common  semiconductors  have  indirect  band-gaps  (the  minimum  in 

c 

the  conduction  band  is  not  over  the  maximum  in  the  valence  band) .' 
To  excite  states  in  these  gaps  would  require  phonon  excitation  in 
addition  to  electronic  excitation.  Furthermore,  the  form  of  the 
wave  function  in  the  indirect  gap  is  not  readily  obtainable  from 
a  simple  model. 

Secondly,  the  various  possible  planes  of  a  three-dimensional 
surface  can  lead  to  surface  states  between  some  bulk  bands  and 
not  others.  The  surface  charge  will  also  be  delocalized  over  the 
surface  plane,  in  contrast  to  only  on  the  other  end  atoms  of  a 
linear  chain. 

Finally,  the  effect  of  such  phenomenon  as  surface  relaxa¬ 
tion  and  reconstruction  are  not  included  in  our  one-dimensional 
model.  Surface  states  can  also  be  modified  by  the  presence  of 
faults  and  adspecies. 

Some  of  these  problems,  including  effects  of  higher  dimen¬ 
sions  and  the  existence  of  adspecies,  are  the  subject  of  con¬ 
tinuing  research.  Dynamical  processes  such  as  laser-stimulated 
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desorption  and  adsorption  are  under  consideration. 
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Figure  1.  Dispersion  relationship  in  complex  crystal  momentum 
space  (k+itc)  for  a  finite  linear  chain.  The  valence,  surface 
and  conduction  bands  are  labeled  V,  S  and  C,  respectively. 


Figure  2.  Absorption  cross-section  for  surface  states,  o,  in 
2 

A  versus  the  frequency  of  the  exciting  laser  radiation. 
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